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1 Introduction 



The nested Bethe Ansatz (NBA) is probably the most sophisticated algebraic construction 
of eigenvectors for integrable lattice models. It appears in models where the underlying 
quantum group is of rank larger than one. 

In the context of the algebraic Bethe Ansatz ]]J the NBA for the A n _i trigonomet- 
ric and hyperbolic vertex model is given in |2j where eigenvectors and eigenvalues are 
obtained for periodic boundary conditions (PBC). In ref. || the NBA for the Sp(2N) 
symmetric vertex model is given and in ref.|3[] the NBA for 0(2N) symmetric vertex 
model is constructed (always with PBC). Although the NBA equations has been pro- 
posed for all Lie algebras H for PBC, no general construction is yet available for the 
corresponding eigenvectors. 

For fixed boundary conditions, the algebraic Bethe Ansatz is known for the six vertex 
model and for the susy t-J model f7j. 

We present in this article the NBA construction of eigenvectors and eigenvalues for the 
A n _i trigonometric and hyperbolic vertex model transfer matrix in the fundamental rep- 
resentation with fixed (SU q (n) invariant) boundary conditions (b. a). That is, boundary 
conditions determined by matrices K ± which satisfy the integrability condition together 
with R{9) 

The NBA is necessary to solve vertex models associated to Lie Algebras with rank 
n — 1 > 1. For the six- vertex (Ai) model, the algebraic Bethe Ansatz gives the transfer 
matrix eigenvectors as products of creation operators of pseudoparticles B{6) acting on 
the ferromagnetic ground state. When the rank n — 1 of the associated to Lie Algebra 
is n — 1 > 1, one finds more than one creation operator for pseudoparticles : B a (6), [2 < 
a < n]. Hence, as Bethe Ansatz for the transfer matrix eigenvectors, a linear combination 
of £> a 's acting on a ferromagnetic ground state state and summed over the indices a is 
proposed. Then one should find the coefficients in such linear combination from the 
eigenvalue condition. Surprisingly enough, these coefficients turn to obey an eigenvector 
problem analogous to the original one but with a new transfer matrix. This new transfer 
matrix is built from statistical weights obtained from the original ones deleting the first 
row and column. This procedure can be iterated as many times as necessary till one 
arrives to a one-by-one transfer matrix. Then the problem is solved in the sense that 
reduces to a set of algebraic equations : the nested Bethe Ansatz equations (NBAE). 

The use of fixed boundary conditions seriously complicates the resolution task. First, 
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the commutation relations between the pseudoparticle operators B a {9) and the transfer 
matrix are much more involved than for PBC and generate therefore many more terms 
when the transfer matrix is applied on the NBA vectors. Second, the structure of the 
unwanted terms generated then is much richer. There appear new algebraic identities 
that were trival in the periodic case and have to be proved now [see eq. (flTf)]. 

In sec. 2 we review the A n _\ trigonometric and hyperbolic vertex model with fixed 
b. c. and its associated SU q (n) invariant spin chain. In sec. 3 we present the NBA 
construction of eigenvectors for this model when fixed boundary consitions are chosen 
and derive the NBAE. This purely algebraic construction is valid for lattices of arbitrary 
size N. We try to keep our presentation as pedagogical as possible. Some calculations are 
given in the Appendices. In sec. 4 we solve the NBAE in the thermodynamic limit. We 
explicitly find the contribution to the free energy of the A n _i trigonometric and hyperbolic 
vertex model due to the presence of the boundaries. From it, we derive the boundary 
contribution to the ground state energy of the SU q (n) invariant hamiltonian. 

2 Construction of the SU q {n) invariant spin chain 

The nonzero elements of the A n _i i?(#)-matrix in the fundamental representation can be 
written for the ferromagnetic regime as: 

RfM = 

R&(0) = 
R2(8) = 

1 < a, b < n 

All other elements are zero. For n = 2, eq.([j]) reduces to the six vertex R-matrix up 
to a gauge transformation f2|. 

The weights in the antiferromagnetic and gapfull regime follow from eq.(|I|) upon re- 
placing 7 — > —7 + Z7T. 

In the gapless and antiferromagnetic regime the R-matrix takes the form: 

T>ab/n\ _ S1I1 7 idsignja-b) / i . 

Kab{&) ~ sin( 7 -0) e ' ' 



sinh7 



sinh(6 l + 7) 

sinh 9 
sinh(# + 7) 



) sign(a-b) 



a 7^ b 



2 



sin 9 
sin(7 — 9) 



(2) 



1 < a, b < n 



Gapfull and gapless antiferromagnetic regimes are related by the transformation: 7 
ry, 9 — ► 

In all regimes, -R(#) fullfils the Yang-Baxter equation: 



[l®R{9-9')}{R{9)®l]{l®R{9')} 
[R{9')®1][1®R{9)][R{9-9')®1\ 



(3) 



The R -matrix (|T]) does not enjoy P and T symmetry but just PT invariance. It is 



not crossing invariant either but it obeys the weaker property || [12 

,1 -1 



S l2 (9) 



n *2 



L(9,j)M 2 S 12 (9 + 2ri)M^ 



where S = PR (Pj£ = Sl5 J k ) and r], L, M are given by |T 



(4) 



7/ = 

M ab = 
L(9,j) - 

Also this R-matrix obeys: 



n 
2 7 



5 ab e {n - 2a+lh l<a,b<n 
sinh(6 l + 7) sinh[# + (n — 1)7] 
sinh(6 l ) sinh(6 l + 717) 



(5) 



R(9)R(-0) = 1 (6) 
We will consider in this paper boundary conditions defined by the K-matrices : 



KM 



>-2a+ih sinh(2g + 7) 



sinh(26 l + 717) 



(7) 



1 < a, 6 < n 



3 



for the right and left boundaries, respectively. They are solutions of the equations 
10: 



R(9 - 6')[K-(6) ® 1}R{9 + 9')[K-{9') ® 1] 

[K~ (9') <g> 1}R(6 + 9') [K~ (6) ® l]i2(0 - 0') (9) 



i2(0 - / )lfi t -(0 / ) tl M 1 - 1 i2(-0 - 9' - 2r ] )K+(9) tl M 2 
= K+{9) tl M 2 R(-9 -9' - 2 v )M- 1 K+{9') tl R{9 - 9') (10) 

Notice that the solutions to these equations can be multiplied by an arbitary function of 
9. These functions were chosen in equations (0), (§) in order to have the term proportional 
to A in the transfer matrix with coefficient equal to 1 (see eq . (f26|) ) . 
The Yang-Baxter operators T a b(9,Co) are defined as usual: 

T ab (9,u)= Yl t aib (e + w N )®t aaax {e + w N - 1 )®...®t aaif _ 1 (e + wi) (11) 

ai,...,ajv_i 

where N is the number of sites, to = (un, Wjv-i, • • • , k>i) and Ui (1 < i < N) are 
arbitrary inhomogeneities. These operators obey the relation: 

R{9 - 9') [T(0, Co) <g> T{9', Co)} = [T(9 f , Co) <g> T(9, Co)}R{9 - 9') (12) 
The row to row transfer matrix for periodic boundary conditions is given by: 

r(9,Co) = Y J T aa {9,Co) (13) 

a 

For fixed boundary conditions described by the matrices K ± (9), one uses the Yang- 
Baxter operators U a b(9,Co) defined by|J: 

U ab (9,Co) =Y,T ac (9,Co)K; d (9)T^(-9,Co) (14) 

cd 

Here T^ l 1 (9, Co) is the inverse in both the horizontal and vertical spaces. That is: 

J2T ab (9,Co)T bc 1 (9,Co) = 16 ac (15) 

6 
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Where 1 is the identity in the vertical space. 
The YB operators U a b(6,Co) fulfil the YB algebra: 



R{6 - 6') [17(0, Co) ® 1}R{6 + 9') [U{6', Co) <g> 1] 
= [U(6\ Co) ® l]i?(6> + 6') [U(G, Co) ® l]fl(0 - 0') (16) 

The fixed boundary condition transfer matrix is defined as: 

t(e,Cu)=J2KMUat,(0,Co) (17) 

ab 

where Co = (uon,uon-i, ■ ■ ■ (see figure A). Thanks to eqs. (f|), (p~0|) , fll6|) and (p~7|) the 
t(0, (D) form a one parameter commuting family: 

[t(9,uo),t(6',oo)] = (18) 
Furthermore, these transfer matrices built with K given by (0)-(@) commute with 



the SU q {n) generators as shown in refs. [IT], p^j . 



The SU q (n) invariant hamiltonian associated to this transfer matrix is given by JTTJ: 



H = 

N—l n r—l s s r— 1 

E{ e (n n ( j n 0+1) + n ( j n ij) n + 

j=l ., (=s £=) — 1 Z=r— 1 Z=s 

r, s = 1 

r > s 

£ s ( n _ r )(h^h^ + h^h^) + ^r(n- r)/W +1 >] + 

n . r=l 

r, s = 1 
r > s 

£ s(r -s)(n- r)(h^h^ +l) ~ hfh^)} + 
r,s = 1 
r > s 

s ^J2r(n-r)(hW-hV) (19) 

where we have ommited a term proportional to the unity operator. Here N is the num- 
ber of sites, J ; + = eu + i, Jf = ei+u and hi = en — ei + u + x are the SU(n) generators in the 
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fundamental representation with (ej m )y = 8u8 m j. It is easily seen that this hamiltonian 
coincides for n = 2 with the SU q (2) invariant one, discussed in [[J-[|10 . 



3 Nested Bet he Ansatz for the open SU q (n) invariant 
transfer matrix 

In this section we give the NBA construction for the A n _i vertex model with open bound- 
ary conditions. 

To make contact with the known case n = 2 is convenient to work with slightly modified 
local vertices: 

[t ab (9)U = R%(e - 7 /2) (20) 
It is also convenient to introduce the notation: 



A{6) = U n (0) 

B a {&) = U la (9) (21) 
V ab (6) = U ab {6) 

2<a,b<n 

The Yang-Baxter algebra fulfilled by these operators follows by inserting eqs. (|l|) and 
(PI) in eq. ( |16|) (see appendix A). 

Actually it is more convenient to work with the operators: 



Ad{d) = — tt^ -[e 2e ^ S mh(29)V bd (e)-smh 7 5 bd A(9)] (22) 

smh(20 — 7) 

The commutation relations are then given by: 

smh(fl + V ) smh{V — o ) 



sinh(2#) sinh(# - 6') 
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sinh 7 e^'sinh(2^- 7 ) 4w ^ c( ^ 



Ad(d) B c {9') 



sinh(20) sinh(0 + 0' 
sinh(0 + 0' + 7) sinh(0 - 9' + 7) 



sinh(0 + 0') sinh(0 - 0') 
i?( 2 )(0 + 9')% BF>{0 - 6')f d B e {9') V gi {9) (25) 

sinh7e e - e, sinh(20 + 7) R (2) (n n )g b m) Am 
+ sinh(0 + 0') sinh(20) * V 9 U*i>WW) 



sinh 7e e e ' sinh (20 + 7) 
sinh(0 - 0') sinh(20) 



Where i?^(0)^ is the original R matrix but with indices 2 < i,j,k,l < n. The 
first term of the last equation may be seen as the building block of a transfer matrix 
of a problem with n — 1 states per link, inhomogeneities 9' and local weights given by 
[t ( 'ab(9)]cd = R^(9)^ with indices going from 2 to n (notice the change 0^0 + 7/2 in 
the local wheights from eq. (|20D ). These commutation relations reduce to those of [JTOj for 
the case n = 2 after cancelling the exponentials by an appropiate gauge transformation. 
Our aim is to build eigenvectors of the transfer matrix t(9,u) (defined by eq. (|TT|) ) . We 
find using eqs. (O) and (f2~2|): 



t(9,Cu) = Y, K a b (0-l/2)U ab (9,Cu) 



06=1 



smh(20 + (n — IJ7J ^ 



Where K^ + {9) is obtained from eq.(^) by making n — > n — 1 and reordering indices 
such that they run from 2 to n. K^ 2 > + {9) is the i^ + -matrix for the reduced problem with 
local vertices [t^(9)] cd = R^(9) b c i 

It is easy to find an eigenstate of t(9, u), the so called reference state || 1 > given by: 
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1 >= 



/ 1 \ 





(27) 



This ferromagnetic state is an eigenvector of both A(8) and T>dd{9) (2 < d < n) with 
eigenvalues: 



A(0) || 1 > = || 1 > 
V dd {6) || 1 > = A_(0) || 1 > (28) 

Where, see [appendix B]: 

a m= 2e yr smh(g + oy- 7 /2) sinh(g - ay- 7 /2) 
" 1 J sinh(# + w< + 7/2) sinh(# - uo t + 7 /2) 1 ' 

In addition, we find [see appendix B] that: 



%{9) || 1>=0, i^j 

U al (6) || 1 >= 0, a > 2 (30) 

Hence, only the a (#)'s acting on || 1 > give some nonzero vector, not proportional to 
the || 1 > itself. 

Therefore, in order to build generic eigenvectors we repeatedly apply operators Bj(/ij) on 
the reference state || 1 > and consider linear combinations. That is : 



*= J2 ^ 1 ---^4(Mi)--A(/v)||1> = &(m)®...®&MX\\l> (31) 

2<ij <n 

Here \i\ . . . fi r and are arbitrary numbers. They will be constrained by requiring 

^ to be an eigenvector of t(9, Cj). We can assume \& to be 6 independent thanks to eq. (|i8|) . 
Our strategy goes as follows. Since t(9,£u) is a linear combination of A(9) and V a a{9) 
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(see eq.(^)) we can apply separately each operator to Then we will use the commu- 
tation rules (^) and (^) to push the operators A{9) and V aa {9) through the BiAfij) till 
A and T) ah reach || 1 >. We use then eqs.(^) and (^). Many terms arise in this way. 
they can be classified in two types: wanted and unwanted. 
Wanted terms are those containing the original vectors: 

BM...BM || 1 > (32) 
Unwanted terms are those where some argument fij is replaced by 9. That is, terms 
arising from the second and third terms in the eqs. ([24]) and (p5|) . These terms are called 
"unwanted" since they can never be proportional to \l/ (here 9 is an arbitrary complex 
number). 

The wanted term in A(9)ty easily follows by repeatedly using the first term in eq. (|24|) . 
We have: 

wanted term in A(9)9 = \[ ^+^-7)^-^-7) ^ 

v ; fJi smh(0 + fjij) smh(0 - //,) v ' 

An unwanted term where B(9) replaces B([ii) follows by using the second term in the 
rhs of (|24j ) when commuting A{9)Bi 1 {^ii) and the first term in (^) for the subsequent 
commutations A(9)B i:j (fij) (2 < j < r). We find: 

sinh7sinh(26 l — 7)e 6l ~ An Jy sinh(/zi + fij — 7) sinh(/ii — fij — 7) 
sinh(26') sinh(# — /ix) j-J~ 2 sinh( / u 1 + fij) sinh( / u 1 — fij) 

B i MBM---&irMX i ^ || 1> (34) 

This calculation was rather simple because B^ (pii) was the first operator from the left. 
Now, we can find the other unwanted terms by pushing the respective i3's to the left using 
the following cyclic symetry implied by eq. (|101|) in Appendix A: 



B(ni)®B(n 2 ) <8> • • • <8> B(n r ) = 
B{n 2 ) ® £(/i 3 ) ® . . . ® B(u r ) ® B(iii) r {2) (Mi, /i) (35) 



where: 



a=2 
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and Tf$(9,jl) is given by eq.flll) with r sites, indices eij running from 2 to n and local 
weights [t ab (0)]ij = Rfo(8). That is, T^ b (9, jl) is the Yang-Baxter operator for a restricted 
model with periodic boundary conditions, one less state per link than in the original model 
and inhomogeneities fx = (ji r , . . . (see figure B). Notice that the inhomogeneities in 
this restricted model are given by the parameters fij of the BA vectors (pTj) . 
From now on, as in formula (^), indices corresponding to lines carrying identical inho- 
mogeneities will be contracted, (see for example figure D). 

Equation (^) tells us that the cyclic permutations jii — > /ij+i followed by the action of 
r^(n,p) leaves \1/ invariant. This property obviously generalizes as : 



B(m) ® . . . B(n k ) ® . . . ® B(nr) = 
&(li k ) ® B(v k +i) ® • • • ® B(Mfc-i) r£\ . . . n (2) (37) 

where rj 2 ' 1 = T^ 2 \^j,fl). 

Using this we can predict the form of the general unwanted term where B{9) replaces 
B{nk) by looking at eqs. (g§ and (|37|). We find: 



sinh 7 sinh(26» - 7) e e ~^ k e e ~^ 



sinh(26») ~ sinh(6» - /i fc ) 

t~t sinh(/i fc + /ij - 7) sinh(/i fc - - 7) 
sinh(/x fe + ^) sinh(// fe - jjj) 

J = 1 

B{6) ® B(fx k+1 ) ® . . . <g> S(// r ) ® ® . . . ® B(^ fe _i) (3? 

rS...^X || 1> 



The third term in the right hand side of eq.fl24|) produces another kind of unwanted 
terms. (Notice that such terms are absent for periodic boundary conditions f2|). We find 
from the third term in eq.(^) using then the first term in eq. (P5|) (r — 1) times and the 
preceding argument: 



sinh7sinh(26> - 7) ,A e d '^ k 



sinh 26 ^ sinh(0 + 
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yr sinh(/i fc + fij + 7) sinh(/i A . - fij + 7) 
sinh(/x fc + Hj) sinh(/x fc - fij) 

J = 1 

B{0) ® B(fi k+1 ) ® . . . ® B(^) ® ® . . . ® B(/i fc _i) (39) 

t^ k ;p)r^...r^ || 1 > X 

where t^(fik', p) is a transfer matrix like in eq.([Tj]) but for a reduced model with n — 1 
states per link, indices running from 2 to n, local weights given by [t^ (#)] c d = R^H^ct 
and inhomogeneities /2 = . . . , fix, fi r , . . . , //fc+i, //&). We have also used (see figure 

C) that: 

n 

£fl (2) (20)5*2 )+ (0) = <k (40) 

and the cyclic symmetry argument above (eq. (|37|) ). Notice that this term is absent in the 
periodic case 0. 

This completes the analysis of A(6)^. 

Let us now compute the action of (see eq. ([26|) ): 

on 

As before, wanted and unwanted terms appear. The wanted term follows by using repeat- 
edly the first term in the right hand side of eq.(|2"5p when T> aa (9) is commuted through 
the Bi/ij). We find: 

wanted term in ~ 7 | e" 29 £ ^{^{0) * = 

sinh(20 + 7) ^ 

e - 2 e sinh(2g - 7) X sinh(0 + ^ + 7) sinh(fl - ^ + 7) 

sinh(2# + 7 ) /J sinh(0 + fij) sinh(0 - fij) ~ { ' 1 1 

BjM...BjM \\i>t^\drtit± xil --- ir 

with fx = (fi r , fl k+1 , fl k , . . . , fix). 

We have collected the R-matrices from the first term in the right hand side of fl25|) into 
t^ 2 \6\p) as in the case of eq.fl3D|) (notice the change 6 — 7/2 — ► 9 with respect to the 
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original problem (20)). 

We see that this wanted term will be proportional to \& if the coefficients X n " 4r form an 
eigenvector of the reduced transfer matrix t^ 2 '{9\ p). That is, if we require: 



t {2) (6;jl)X = A (2 \6;j2)X (43) 

The unwanted term coming from the second summand in eq. (p^| ) follows by the usual 
symetry argument after using eq. (|40|) and the first term of eq.(^4|) r — 1 times. This gives: 



sinh 7 sinh(2fl - 7) ^ e' 6 '^ 



n 

j = i 



sinh(26') ^ sinh(# + fi k ) 
sinh(/x fc + fj,j - 7) sinhpfc - - 7) 
sinh(yu fc + fij) sinh(/i fe - fij) 



B{6) ® B{ii k+1 ) ® . . . ® B(fx r ) ® ® ® . . . ® B(/i fc _0 (44) 

ra...r< 2 »X||l> 

Notice again that this term is absent in the periodic it happens with eq. (j39|) . 

The last term coming from the action of eq. (|4"ID in \l/ follows using the third and first 
terms of eq.(p5|) and the identity (fiTf). We get : 

sinh7sinh(26 l — 7) e ~ e -Hk 

2^ • ^(n 



sinh 2$ ^ sinh(0 — 
sinh(/i fc + fij + 7) sinh(/i fc - fij + 7) 
sinh(/i fc + fij) sinh(/i fe - fij) 



n 

j = i 

® ® . . . ® ® ® . . . ® B(^fc_i) (45) 

with /i = (/ifc-i, . . . , /ii, /i r , . . . , /ifc+i, /i fc ). 
The sum of wanted terms reads from eqs. (|33|), fl42|) and (1501): 



wanted term in t(6, lj)^ = 
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r -Jj sinh(6 l + fij — 7) sinh(0 — fij — 7) 
^ sinh(6» + sinh(6 l - /^) 

sinh(2fl - 7) SL sinh(g + ^ - 7 /2) sinh(fl - ^ - 7 /2) 

sinh(2# + 7) W sinh(# + ^ + 7 /2) sinh(fl - ^ + 7 /2) [ ' 

fj sinh(g + fXj + 7) sinh(fl - //j + 7) (2) 

y sinh(0 + // j )sinh(0-// i ) ^ ,/iJJ * 



where we have also used eq. (|29|) . The term in brackets gives us the eigenvalue of the 
initial problem in terms of A^ 2 \6;fi), the eigenvalue of the reduced problem defined by 
t (2) (6, ft) with n - 1 states per link and local weights [^(9)]^ = R® bd (9). 
Before summing the unwanted terms we use the identity: 

*< 2 > fa, H)tI% . . . rf > = r& . . . r?¥V fa, ft (47) 

where ft = fa-i, ■ ■ ■ , A*i> Vr, ■ ■ ■ , Vk+i, Hk) and fi = fa, fi k+1 , fi k , fi k _ x , . . . , 
This identity tells us how to move cyclically the inhomogeneities in the open chain. Al- 
though this is a trivial rotation in the periodic case, this is not the case for open boundary 
conditions. The proof for three sites is in figure E, it uses the Yang-Baxter eq.(D), the 
property fllCf) and eq.(f|). Notice that it is enough to prove that: 

t^fa,M 2) =A 2) t^fa,fx) (48) 

with fl = (//!, fi r , . . . , /i 3 , /i 2 ) and fi = fa, . . . , fi 3 , fi^, fix), (see figure D). The proof for 
an arbitrary number of sites is straightforward using repeatedly what was used for three 
sites. 



Using this property one obtains from eqs. (pq), (pM), (14) and (EH): 



sinh (2$ — 7) sinh 7 >J — 

r 

[ n 



1 



k J t smh(6 + fi k ) sinh (6 - fi k ) 
sinh(/i fc + fij - 7) sinh(/i fc - fij - 7) 



J = l 

r 

n 

j = i 



sinh( / u fc + fij) sinh( / u fc - fij) 



smb. fa + \ij + 7) sinh(/i fc - fij + 7) 
sinh(/i fc + fij) sinh(/i fc - fij) 
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i\ sinh(/i fc + Uj- 7/2) sinh(/i fc - ay- 7/2) (2) 

ii sinh(/i fc + Wi + 7/2) smhOi* - w, + 7/2) 1 } 

B{9) ® ... ® 6(11^)4% . . . r[ 2) X (I 1 > 



where we have used 
The unwanted terms have to be zero if \I> is to be an eigenvector. 
In summary, we find the two following conditions to be satisfied : 



t^(9;jl)X = A®(9;ji)X (50) 

A ( 2 ) / ~ x A sinh(/i fc + Ui + 7/ 2) sinh(/i fc - u t + 7/ 2) 
H[ smh(/i fe + LUi - 7/2) smh(/i fc - u)i - 7/2) 
-pr sinh(/i fc + fj,j - 7) sinh(/i fc - fij - 7) 
sinh(// fc + fij + 7) sinh(// fc - fij + 7) 

J = 1 

It is easy to see that eq. ([[I]) ensures the analyticity of the wanted term flilf ) as a 
function of 6 1 for 9 = 1 < j < r. 

We have reduced the original problem of N sites, n states per link and local weights given 
by (|2~Cf ) to a problem of r sites n-1 states per link and local weights [t a b(9)] cd = R^ 2 \9) b c ^ 
with inhomogeneities \i\ . . . \x r . 

By analogy, we propose the following ansatz for the coeficients X^ 1 ' = X %1 " ,tr : 

x« = x® ® . . .®m^fi^) 11 iw > (52) 

where || l^ 2 -* >= ®fLi II 1 > an d || 1 >^ is a n — 1 component vector with the first 
component equal to one and the rest vanishing, and (if^ = 1 < i < r = p\. 
This argument can be repeated as many times as necesary till the dimension of the 
vertical spaces reduce to one. We get in this way a sequence of Bethe Ansatz, each of 
them contained in the previous one. That is, a nested structure emerges. 

It is important to remark that the spectral parameter and the roots of the Bethe 
ansatz suffer in the course of the construction a change 9 — > # + 7/2 from a Bethe Ansatz 
at a given level to the next, and that the roots at each level are the inhomogeneities for 
the next level. This can be seen looking to the first term of the commutation relations 
( p.00p in Appendix A. Then one obtains: 
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*<* +1 >(0,/2 (fc) )* (fc) = A (w) (^ w )I (t) 



A( *), fl -(fc-m _ ft Sinh[d + ^ + (fc ~ 2)7] Sinh(g ~ ^ ~ 

^ ' W -urn , ( fc ) , n i • i (a ( k ) 

j=i sinh|0 + fij + (k — IJ7J smh(0 — //j ' 

sinh[2# 4- (jfe - 2W1 



7) 



4 

Pfc-i 

r 

3 



anh[20 + (k — 2)7] 
sinh(20 + fry) 
n sinh[fl + /if ~ 1} + (fc - 2)7] sinh(fl - /if ~ 1} ) 
M sinh[# + /if " 1} + (fc - 1)7] sinh(# - /if ~ 1} + 7) 
fr sinh(g + /if + fry) sinh(fl - /if + 7) ^ & 
ji sinh[fl + /if + (k — 1)7] smh(9 - /if) 
1 < jfe < n - 1 , /if = ujj + 7/2 , A (n) (#, /x^" 1 )) = 1 

with /if ^ obeying: 



A( fe+1 )(/if ) ,/i«) = 

^ sinh[/if } + /if - x) + (fc - 1)7] sinh(/if } - /if ~ 1} + 7) 
M sinh[/if } + /if ~ 1} + (jfe - 2)7] sinh(/4 fc) - //f 
|| sinh[/if } + /if ) + (k - 2)7] sinh(/if ) - /if } - 7) 
sinh(/if + */f + fry) sinh(/xf — uf^ + 7) 

J = 1 

J 7^ « 

Using the recurrence formula ([54]) we find for the eigenvalue of £>) : 



ACOrfl „(0)x _ fr siuh(g + /if - 7 )sinh(g-/if) 

iV (tf,/i J — 11 . {0) . {0) 

j=i smn(0 + //] ; ) smh[V - fiy + 7) 
" sinh(2fl - 7) sinh(2fl) 

V"- 1 smh[6 + /if _1) + (jfe - 1)7] sinh(# - /if _1) + 7) 
sinh[# + /if ~ 1} + (jfe - 2)7] sinh(# - /if -1 *) 



^ sinh[2# + (jfe - 2)7] sinh[2# + (jfe - 1)7] 



i=i 

^ sinh[# + /if } + (jfe - 2)7] sinh(# - /if } - 7) 
II 



( fe ) , /-fc iu„;„U/i ,.( k h 



1 sinh[# + /if + (jfe - 1)7] sinh(fl - //} 



15 



Where in the last term the product over p n is substituted by 1. 

Let us now derive the Bethe Ansatz equations for the parameters /4 (1 < i < Pk , 1 < 
k < n- 1). 

Changing by + 1 in eq. (0) and setting 9 = /4 yields: 



A (fe+1 Vf\/i (fc) ) = 

^ sin%fj + /if +1) + (fc — 1) 7 ] sinh(/.f } - £ +1) - 7 ) 
3 -=i sinh(/4 + /if +1 ^ + fcy) sinh(/4 — {if ) 
since the second term of eq.(|54|) vanishes for this value of 9. Equating now eqs. (|55|) 
and (151 



j^j- sin(i/$ + z^- + ry) sin(i^- — i^- + 27) 
sin(V} + z/] fc ^ — 27) sin(i/j — i/^ — vy) 



3=1 



sin 






+ i-y/2] 


sin( 


(fe) 
f i 


-z/f- 1) +z 7 /2) 


sin| 


- (fe) 


1 (fc-1) 


- i-y/2] 


sin( 


(fc) 
f .■ 


- vf^ - .7/2) 


sin 


- 0) 


+ +1) 


+ 27/2] 


sin( 


(k) 


-*? +1) +,7/2) 


sin 


- (k) 


+ vf +1) 


- 17/2] 


sin( 


(k) 


- ~ *7/2) 








1 < 


i < 


Pk 


, 1 < jfe < n- 1 



(58) 



where we have set: 



= iuf } - (ife - 1)7/2, 1 < jfe < n - 1, 1 < j < p k (59) 
The function A^\9, fx^) must not be singular at the points 9 = jjL k - (1 < j < p^ , 1 < 
k < n — 1) since the finite dimensional matrix t(8,u) is an analytic function of 9. One 
can see that if the previous Nested Bethe Ansatz Equations (|58"D are satisfied by the fij 
the residue of A^(9, li^) at 9 = Lip and 9 = —j^P — (k — 1)7 vanish. 
The NBAE for the gapless regime follow from eqs. fl5*B"|) by replacing 7 — > — 27, z/j^ — > zz^. 



4 Analysis of the Bethe Ansatz equations 

In this section we investigate the solution of the NBAE associated to the quantum group 
covariant NBA with the inhomogeneities at the first level fixed to zero for simplicity. 
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We will relate these equations with those of the periodic case by means of the change of 
variables, (see fT0|j ): 



%P k - 8+ i = ~"i k) (60) 
1 < s < pk ; 1 < k < n — 1 



In this way the NBAE can be written for the gapless case as: 



2 ^ sinh(A{ fc) -Af +i 7 ) 
M sinh(Aj fc) - Af - vy) " 

sinh(Af } + zt/2) sinh(Aj fc) - »(tt - 7 )/2) (gl) 

sinh(Af } - z 7 /2) sinh(Af } + z(tt - 7 )/2) 

2 ^ sinh(A ; (fc) - Af -1) + ry/2) 2p J^ sinh(Aj fc) - \f +1) + z 7 /2) 

M sinh(Af > - Af - 1] - vy/2) M sinh(Af } - Af +1) - ry/2) 
1 < k < n - 1; 1<1 <2p k 



[In the gapless regime, the statistical weights are given by eq.(|2]) ]. 
These equations are like the NBAE for periodic boundary conditions on 2N sites with an 
additional source factor, (see 0). In addition we have the following constraints on the 
roots A^: 

(i) the total number of roots is even {2p k ) at every stage and are symmetrically distributed 
with respect to the origin according to eq. (|60|). 

(ii) there is no root at the origin A^ = at every stage due to the fact that Bj k \9 = 
—(k - 1)7/2) = (see eq. @). 

As usual, we take logarithms in eq. (Wfh , yielding: 



2 Pfc+i 2p fc 

E ^-A^.T/^-E^-^T) 

3=1 3=1 

+ E - Af-^,7/2) + $(Af } , 7 /2) - <S>(\?\ (tt - 7 )/2) (62) 

3=1 

= 2vr/f } 1 < jfe < n - 1; 1 < / < 2p fc 
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Where 7 Z are integers and 



${z,i) = i log 



sinh(i7 + z) 



(63) 



sinh(i7 — z) 

We will consider the thermodinamic limit of eq. (|6TD . One can now introduce a density 
of roots at every NBA level: 



p«(Af)= lim — !- , 

N-,00 AT \W _ \('h 



mu - A )) 



Defining now the counting functions as: 



(64) 



2PM 



(fc+1) 



2vriV 



2p fe - 



i=i 



+ £ $(A-Af- 1 \ 7 /2) + $(A )7 /2)-$(A,( 7 r- 7 )/2)] 
i=i 

1< Jfe < n- 1 



(65) 



And using that: 



N 



(k) 



j{k) _ j(k) _ 1 r 



h=l 



33h(k) 



(k) 

Where is the number of holes at level k. One can see for N — > 00 that: 



(66) 



r(k) 



(*)/ 



rfA 



Aj + l - Aj 



(A i+1 -AjOiV ~ iV AT 



(67) 



Where the term is produced by the hole at A = at every NBA level. In the limit 
of large N we have: 



i im >T7E/(Af)= / d\f(\)pM(\) 

N^oo iV J~ [ J J-00 



(6? 
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Taking the derivative of eq.(^) with respect to A and using eq. ([67]) , we obtain integral 
equations for a^ k \X). Let us start with the antiferromagnetic ground state. That is, no 
holes besides A = 0, and no complex solutions. We find: 



^ 1 />oo 

a k (X)-J2 / d^K km (X - //)cr m (//) 

m=l- y -°° 



' $'(A, 7 /2)--l-$'(A,(7r-7)/2) (69) 



2vriV v " ' 2nN 

$'(A,7/2)-^E^(A) 

(0) 



Where we set all the inhomogeneities equal to zero at the first level ( LUi = — > A 
0). The kernel K km (X) reads: 

2vK km (X) = $ , (A, 7 /2)(5 fcjm+1 + 5 fejm _ 1 ) - $'(A, 7 )4 m (70) 

This linear integral equation can be solved by means of the resolvent R mn (X) given by 
the solution to the equation: 

n— 1 

Rik(T - X)[8 km 5(X- /j,) - K km (X- /i)]dX = 5(t - ii)5 lm (71) 

fc=i J -°° 

It is convenient to Fourier transform these quantities: 

—e* kX R rnn (k) (72) 

-oo Z7T 

^(A) = /^§^Mfc) (73) 
The solution to eq. ([fl]) is then given by 0: 

Rn' (2a;) = sinh ( 7ra; ) sinh[7x(re - />)] sinh^q^) 
sinh[x(7r — 7)] sinh(7xn) sinh(7x) 

where /> = max (I, I') and /< = min(l,V). We then obtain for the derivative of the 
counting functions: 
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, 2sinh[ 7 fc(n-/)/2] 
K) sinh( 7 fcn/2) 



i _ i ( ^W^.-7M ) g ^ (75) 

iV iv smh(/c7r/2j 



One can see that this result reduces to the one given in [[TtJ for the case n = 2. 
To compute the physically meanigful quantities only p^>(k) is needed. Using eqs. 
and (|75D is easy to see that: 



. ( = ^MLzJO/g] (76) 

sinh[/c7(n — l)/4] cosh(/c7r/4) 
AT cosh(fc7n/4) sinh[A;(7r — 7)/4] 

We have now the tools to evaluate the free energy of the model in the gapless regime. 
This is given by: 



mi,n) = --J-logA(0) 



N -> oo 



N -> oo 



AT 



-A,, 7 /2) (77) 

~ ^ e _ M sinh[^ - 7 )/2] (1 , 
-oo k smh(/c7r/2) 



Note that we have made 9 — > # + 7/2, that is we have returned to the local weights 
where i(0, tD) oc 1. 

Using the expresion for p^ l \k) (eq. (IT7|)) in ( |?8| ) the final result for the free energy is: 

f(0, 7, n) = 4 f ^ sinh(2xfl) ^MlLZ DJ sinh[*7(n - 1)] 
jo a; sinh(a;7r) sinr^aryn) 

2 ^ sinh(2x^) COSh ^^ ~ 7)/2 ^ Sinh ^ ~ 1)X7/2 ] 



N Jo x 1 sinh(x7r/2) cosh(xn7/2) 
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The first term here is the known bulk free energy, (see [[J)- The second term is the 
correction produced by the open boundary conditions (that give quantum group invari- 
ance) . 

The ground state energy for the SU q (n) invariant hamiltonian is obtainined by using: 

H = i(0, 0) + (JV - 1) ^ ~ ^ cos 7 (80) 

2 n 

We obtain, (deriving ( f75| ) with respect to #): 



n — 1 , f 00 , sinh[x(7r — 7)] sinhtaryfn — 1)1 , , 

eoo7 = cos7-4sin7/ 1 : . , " . , y \ ^ (81) 

n jo smh(x7r) smh(x7n) 

(n — 1) 2 sin 7 Z" 00 cosh[x(7r — 7)/2] sinh[(n — 1)3:7/2] 

cos I / diOo ~ 



Nn N Jo sinh(x7r/2) cosh(3:7n/2) 

In the special case n — 2, this reduces to the result in fT3|j . 
The surface energy contribution in eq.fl3T|) : 



o. . n—l n . f°° , cosh[x(7r - 7)/2] sinh[(n - 1)3:7/2] 

e 5 7 = cos7 + 2sin 7 / dx 1 . V L ' '! ,1 , (82) 

n Jo sinn(x7r/2) cosh (3:771,/ 2) 

takes a simpler form in the 7 = (isotropic) limit. We find 

c *(o) = + 2 r dx eM ~ x/2] ^ h[( r r 1 W2] , (83) 

n jo cosh(3;n/2) 



This integral can be expressed in terms of elementary functions [15 



e s (0) 



n 



n 



7T 



E cos 

fe=0 



77, 2sin(7r/n) 
2A; + 1 



-ln2 



77 



-71 log 



/2A: + 1 N 

2 — 2 COS 7T 

V n 



} 



54) 



Here stands for integer part of x. 

Let us finally consider the gapfull antiferromagnetic regime. In this case the NBAE 
are solved by expanding in Fourier series, since the NBAE roots are in the interval 
(— tt/2, +7r/2). We write the density of roots as follows: 



00 g2imA 



2tt 



aHm) 
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where ex' (A) obeys a system of integral equations analogous to eq. (|69D : 



where now, 



n-l 



'*~ x r+n/2 

a k (X)-Y, / dnK km {\ - fi)a m (fi) 
m =iV— / 2 

= ^*^.V2) + ^'(A,(„ + oO/2) 

71 iv m=l 



(86) 



$(2,7) = * lo § 



sin(ry + 2) 



sin (27 — z) 

and the kernel iffc m (A) is given by eq.(^). We find as solution of eq . (p6|) : 



a l (m) 



4sinh[7m(n — I)] 
sinh(7mrz) 



+ 



+ 1 ^2 + {[1 + (-1)™] exp(-|m| 7 ) - 1} £ (88) 
Where -R^(m) is the resolvent of eq.(|S6"D in Fourier space. Then, using eq.(^), we 



find 



^ _ 4sinh[7m(n - 1)] + h m 
sinh( 7 mn) iV 



39) 



where 



(— l) m sinh[(n — l) 7 m/2] + exp[— (n — l) 7 |m|/2] sinh(7m) 
cosh [777172/2] sinh( 7 m/2) 
We find upon inserting p l {m) and $(2,7) given by eqs.(^7|)-(|89|) in eq. (|77|) 



(90) 



n 



m=l 



mi , „) = 49(1 - I) + 4 £ e -" ,sinh '-""'"- 1 )' 5inh ' 2 '" 9 ' 

msinh(7mn) 

e~ mi h m sinh[2m0] 

m=l 



+ 



iV 



1)*+ E 



(91) 



The first two terms correspond to the known bulk free energy, (see [0). The second 
term is the correction produced by the open boundary conditions (that give quantum 
group invariance). 
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5 Conclusions 



We have presented the generalization of the Nested Bethe Ansatz to the quantum group 
invariant case. It will be interesting to generalize it to the cases where the K matrices 



are the general diagonal solutions given in [11 



It also remains to study the quantum group properties of the NBA states as the highest 
weight property. Moreover, a rich structure must arise for the reduced models when 7/71 



is a rational number 14 



It would be interesting to study this construction for algebras different to A n -\. That is, 
to generalize the work in refs. ||-|| to open boundary conditions 

6 Appendix A : commutation relations 

We begin putting explicitly all indices in eq. @. This yields (from now on we will supose 
sum over repeated indices): 



MS = R{9 - e>)? f U eg (8)R(0 + 0% f d U, 

n$ = u ae {e>) R(e + e% u fh (6) R{e - efj 



ef Ueg\V)tl{0 + V )fa U hc \V ) = ^ 



As we want to obtain the commutation relations between A(9) = Un(6), T>m{9) = 
U bd {0) and B c (9) = U lc {9) where (6, c, d > 2 ) , we study the equalities M\\ = N£ and 
M$ = N^. This gives: 



sinn (6/ — V') sinh (6/ + (J + 7) 



e-(«-»')sinh(« + «Osi n h 7 



sinh(# - 9') sinh(# + 9' + 7) 

e (6+6 ) ginh^y 



sinh(# + 9' + 7) 



B g (9')V gc (9) (93) 



V bd {9) B c {9') 



sinh(# + 9' + 7) sinh(fl - 9' + 7) 

sinh(# + ¥) sinh(# - 9') 
{R {2 \9 + 9>)f h RP>{6 - 9')f d B e {6>) V gi {9) 

e~^ e ~ e > sinh 7 
~sinh(0-e' + 7) 



R^(9 + 9')f d B g (9) V t 
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e ( e+e ') sinh7 
sinh(6> + 9' + 7) 
e~( e ~ e > sinh7 



.4(0) £ c (0') <U} (94) 



sinh(0 - 0' + 7)' 

Where R^ 2 '(9)^i is the original R matrix but with indices 2 < i, j, k,l < n. We would 
like to have all the i3's to the left of the ^4's in the right hand side of eq.(p4|). For that, 
one substitutes eq. fl9~3| ) in the last two terms of eq. (0). One obtains in this way a long 
expresion that we omit. 

To simplify the calculi of Bethe Ansatz we look now for a linear change of the operators 
such that no term proportional to B g (6')A(8) remains in the commutation of P's and S's. 
The most general linear change would be of the form: 

V hd {9) = aZ(e)V rs (e) + (3 bd (9)A(9) (95) 

With ot r bd (9) an invertible matrix. Plugging this in eq. ( |9~4"D and imposing the cancela- 
tion of terms of the form B g (9')A(9) one obtains: 



<*W) = <*(W5 (96) 
(3(9) /a(9) = - e - 20 sinh7/sinh(2fl + 7) 



We define the operators: 

1 



V bd (9) = _^[e 29 sinh(20 + 7)P M (0)-sinh7^(0)] (97) 

B = ^h(20 + 7 ) 

v ; sinh 20 v ; v ; 



Now, after some work we arrive to: 



A(9) E (9') sinh(g + ff) sinh(0 -ff-j) $ 
A W B * 9 > ~ s inh(0 + 9' + 7) sinh(0 - F) B ° { ° ] A{6) 

sinh. ^yc^~^ sinh. 2$ ■*> 
+ s inh(29 + 7 ) s mh(f)-9') Bc{B) 



sinh 7e »-« sinh 20 SM%(#) (») 



sinh(20 + 7) sinh(0 + 9> + 7) 
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Me)U 0) = ^ ( i^;^ft/; 7) doe) 

sinn(0 + 0' + 7) sinh(0 — 9') 

R {2 \e + e' + 7 )f h R^\e - e't d B e (e') v gi (9) 

sinh7e e-g' s j n h(2g + 2 7 ) (2) 6 ~ ~ 

-sinh(^-^)sinh(2^ + 7 ) * (2£> + 7L ^ Ac( * } 
sinh7e e-e' sin h(2g + 2 7 ) (2) 6 ~ 

+ • r 77 , fl , , — x ■ , foo , — 7 R 'i 20 + Vcd Bg{9) A{9 
smh{0 + 9' + 7) sinh (2 1' + 7) 

Through the transformation : 



9^9-^/2 
9' -> ^-7/2 

Equations (]23j) and @ follow. 
It will be also necessary to derive the commutation relations between the B's. This is 
obtained using the equality Ml} = with (c, d > 2). This gives: 

B c (9') = B g (9') B h (0) R h J(9 - 9') (101) 

7 Appendix B : evaluation of A_(0) for SU(n) 

We will work with a chain of length N (remenber 9 — > — 7/2 for the first level). One 
can easily see that: 



Zu(0) 
2k (0) 



1 > 
1 > 



1 > 



N 

n 



sinh(6» + o>j - 7/2) 
5[ sinh(0 + uji + 7/2) 
1 > ^ 

1> = 0i^ji,j>2 
1 > = 



1 > : = <L(0) || 1 > 



(102) 



25 



Tn(9) 


1 > 


= II 1 > 


f dd {9) | 


1 > 


SL sinh(# 
fjsinh(0 


f u {9) 


1 > 


¥> o 




1 > 




f dl (0) 


1 > 


= 



1 >:=<L(0) || 1 > 



(103) 



We will now evaluate the action of U bd on the reference state. 

A{0) || 1 > = Tu^fnie) || 1 >= 
= T n (d)f n (6) || 1 > 
= II 1 > 



(104) 



U dl (6) || 1> = T dl {9)f n {9) || 1>=0 



(105) 



Where we have made use of eqs. (|102|) and (|103| ). The other elements are not so easy to 
evaluate. One has to make use of the following identity which follows by direct calculation: 



f ld (6) || 1 >= -ei640)T u {-9 - 7 ) || 1 > 
It is also needed to derive from eq.(|T2D that: 



(106) 



T dl (6)T lb (-0- 7 ) = T lb (-e-j)T dl (9) 

e~ 26 sinh 7 



(107) 



+ 



sinh 29 



[T n (-9 - ry)T db (9) - T n (9)T db (—9 - 7)] 



We have in addition: 



U db (9) || l>=T dl (9)T lb (9) || 1> 

We distinguish now two cases: d 7^ b and d = b: 
(i) d ^ b 
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(10* 



U db {6) || 1 > = T dl (6)f lb (6) || 1 > (109) 
cc r dl (0)T 16 (-0- 7 ) || 1>=0 

This can be seen applying both sides of eq.( |107|) to the reference state and using 



eq. 



fin d = b 



U dd {6) || 1 > = T dl (9)f ld (6) || 1 > + 5-(6)5-(6) II 1 > 



[smh(29 - j)5-(9)8-(9) + e- 2e sinh 7 ] || 1 > 



sinh2# 

where we have used eq.( 107|) and the fact that: 



(110) 



To conclude, we have: 



5-(ff)5(-e-i) = 1 



Where: 



A{9) || 1 > 

U dl (6) || 1 > 

V db (6) || 1 > 

B d {6) II 1 > 



1 > 





A_(0)fo || 1 > 



;ii2) 



A-(6) = e 2d 5_(6)5-(6) 
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N 

n 



sinh(6» + uji - 7/2) sinh(6» -u { - 7/2) 



\ smh(0 + uji + 7/2) sinh(6» -uji + 7 /2) 



;il3) 
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8 FIGURE CAPTIONS 



A. Transfer matrix for the open chain with inhomogeneities uj = (ujn, ■ ■ ■ > 

B. Transfer matrix for the periodic chain with inhomogeneities jx = (/in, ■ ■ ■ , A*i)- 

C. Identity fl40|) and its application to construct the open transfer matrix. 

D. Representation of the identity (f48|) . Indices are contracted corresponding with the 
inhomogeneities. 

E. Graphical proof of eq. (|48|) for three sites. We use the following : (a) eq. ( |40D and 
i?(0) = 1. (b) eq. (Q), repeat this step for an arbitrary number of sites, (c) eq. (|J). (d) 
eq. (H), repeat this step for an arbitrary number of sites, (e) eqns. (||) and (flCf). (f) eq. 

©■ (g) R(0) = i. 
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